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HEAT SHIELD CAVITY PARAMETRIC EXPERIMENTAL AEROHEATING FOR A 
PROPOSED MARS SMART LANDER AEROSHELL 


Derek S. Liechty*, Brian R. Hollis f 
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ABSTRACT 

The proposed Mars Smart Lander is to be attached through its aeroshell to the main spacecraft bus, thereby pro- 
ducing cavities in the heat shield. To study the effects these cavities will have on the heating levels experienced by the 
heat shield, an experimental aeroheating investigation was performed at the NASA Langley Research Center in the 
20-Inch Mach 6 Air Tunnel. The effects of Reynolds number, angle-of- attack, and cavity size and location on aero- 
heating levels and distributions were determined and are presented. To aid the discussion on the effects of the cavities, 
laminar, thin-layer Navier-Stokes flow field solutions were post-processed to calculate relevant boundary layer prop- 
erties such as boundary layer height and momentum thickness, edge Mach number, and streamwise pressure gradient. 
It was found that the effect of the cavities varies with angle-of- attack, freestream Reynolds number, and cavity size 
and location. The presence of a cavity raised the downstream heating rates by as much as 325% as a result of bound- 
ary layer transition. 


NOMENCLATURE 

h heat transfer coefficient, h=q/(H aw -H w ), 

(slug/ft 2 / s or kg/m 2 /s) 

H enthalpy (BTU/lbm or J/kg) 

L length of vehicle from nose to base (in. or m) 

M Mach number 

p pressure (lb/in 2 or Pa) 

q surface heat transfer rate (BTU/ft 2 /s or W/m 2 ) 

R model reference radius (in. or m) 

Rjy model base radius (in. or m) 

R n model nose radius (in. or m) 

R s model shoulder radius (in. or m) 

Re unit Reynolds number (1/ft or 1/m) 
r radial distance from symmetry axis (in. or m) 

a distance along surface from cavity center (in. or 

m) 

u velocity magnitude (ft/s or m/s) 

w cavity diameter (in. or m) 

a angle-of-attack (deg) 

S boundary layer height (in or m) 
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(p model orientation angle (deg) 

r] body half angle (deg) 

0 momentum thickness (6^/(pu/p e u e )*(l-u/u e )dy) 

(in or m) 

6 cavity orientation angle (deg) 

p density (slug/ft 3 or kg/m 3 ) 

§ first afterbody cone angle (deg) 

£ second afterbody cone angle (deg) 

Subscripts 

oo freestream static conditions 

aw adiabatic wall conditions 

e boundary layer edge conditions 

FR conditions from Fay-Riddell calculation for a 

hemisphere 

fs full scale dimensions 

w conditions evaluated at the wall 

INTRODUCTION 

The design 1 of the proposed Mars Smart Lander 
(MSL) aeroshell requires that it be attached to the main 
spacecraft bus by way of six structural inserts which 
pass through holes in the forebody heat shield (Fig 1). 
When the lander is separated from the spacecraft bus 
prior to entry into the Martian atmosphere, the bolts are 
severed and retracted, which forms cavities. The pres- 
ence of these cavities in the heat shield during entry may 
result in high, localized heating at the downstream edge 
of the cavities due to flow separation and reattachment 
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Figure 1: MSL cruise configuration. 

within the cavities; and may be accompanied by a 
heating augmentation downstream of the cavities due 
to a change in the state of the boundary layer from 
laminar to transitional or turbulent. The Genesis Sam- 
ple Return Capsule^ (GSRC) also contains forebody 
heat shield cavities, but there are several significant 
differences between the two vehicles. First, the fore- 
body half angle of the GSRC is different than the pro- 
posed MSL. Also, the proposed MSL is being 
designed to enter the Martian atmosphere at an 
angle-of-attack of 16-deg, but most of the data for the 
GSRC collected in Ref. 2 were at 0-deg 
angle-of-attack. 

The goal of the present study is to determine 
experimentally how these cavities affect the heating 
levels on the forebody of a MSL concept. Tests were 
conducted at nominal conditions of Mach 6 in air with 
freestream Reynolds numbers from 2.1xl0 6 to 
7.3xl0 6 per foot. A range of angles-of-attack from 
0-deg to 20-deg was studied. These conditions result 
in boundary layer edge conditions similar to those 
found in flight 3 . Transition correlations have been 
developed that relate the Reynolds number based on 
momentum thickness to the ratio of cavity diameter to 
boundary layer thickness for the wind tunnel condi- 
tions in Ref. 4 and for flight conditions in Ref. 3. 
Three nominal cavity sizes were tested at two radial 
locations. Global heating distributions were measured 
using thermographic phosphors. To complement the 
measurements made in the wind tunnel, laminar, thin 
layer Navier-Stokes flow field solutions generated 
using the Langley Aerothermodynamic Upwind 
Relaxation Algorithm 5,6 (LAURA) were post pro- 


cessed to calculate boundary layer height, momentum 
thickness, edge Mach number, and streamwise pres- 
sure gradients. Only a qualitative discussion of 
boundary layer properties will be included in this 
report. Corresponding computations and correlations 
for this study are presented in the companion paper by 
Hollis 4 . 

EXPERIMENTAL METHODS 

Test Facility 

Aeroheating tests were conducted in the NASA 
Langley 20-Inch Mach 6 Air Tunnel. This is a blow- 
down facility in which heated, dried and filtered air is 
used as the test gas. The tunnel has a two-dimen- 
sional, contoured nozzle that opens into a 20.5-in. by 
20-in. test section. The tunnel is equipped with a bot- 
tom-mounted injection system that can transfer a 
model from the sheltered model box to the tunnel cen- 
terline in less than 0.5 seconds. Run times of up to 15 
minutes are possible in this facility, although for the 
current aeroheating test, run times of only a few min- 
utes were required (models are only exposed to the 
flow for a few seconds). The nominal reservoir condi- 
tions of this facility are stagnation pressures of 30 to 
500 psia (206.8 to 3447.4 kPa) with stagnation tem- 
peratures of 760° to 1000°R (422.2 to 555.5 K), which 
produce perfect gas freestream flows with Mach num- 
bers between 5.8 and 6.1 and Reynolds numbers of 

0.5xl0 6 to 7.3xl0 6 per foot (1.64xl0 6 to 23.95xl0 6 
per meter). A more detailed description of this facility 
is presented in Ref. 7. 

Phosphor Thermography Technique 

Global surface heating distributions were calcu- 
lated using the digital optical measurement method of 
two-color, relative-intensity, phosphor 

thermography 8 ' 11 . Ceramic wind tunnel models are 
coated with a phosphor compound that fluoresces in 
two separate regions (green and red) of the visible 
light spectrum. During a wind tunnel run, the phos- 
phor-coated model is illuminated by ultraviolet (UV) 
light sources, and the resulting fluorescent intensity of 
the model is recorded and digitized through a color 
CCD (charge coupled device) camera. The fluorescent 
intensity is dependent on both the intensity of the inci- 
dent UV light and the local model surface tempera- 
ture. The UV intensity dependence is removed by 
taking the ratio of the green to red intensity images, 
from which surface temperature distributions can be 
determined through prior calibrations. Images are 
acquired before the wind tunnel run and after injec- 
tion of the model to the tunnel centerline during a run. 
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Global heat transfer distributions are then computed 
from these temperature data using one-dimensional, 
constant heat-transfer coefficient conduction theory 1 1 . 

The global phosphor thermography technique is 
now the standard method for aeroheating studies in 
Langley’s hypersonic tunnels used for aerothermody- 
namic studies. The global data obtained using this 
method can be used to identify the surface heating 
effects of complex three-dimensional flow phenom- 
ena such as transition fronts, vortex structures, and 
shock interactions which are difficult to examine 
using conventional discrete-sensor methods such as 
thin-film resistance gages or coaxial surface thermo- 
couples. 

Test Model Description 

In order to manufacture ceramic test models, 
rapid-prototype, stereolithographic (SLA) resin mod- 
els were first fabricated based on surface geometry 
definitions in electronic data files. Wax molds of the 
resin models were made, and then a patented 12 silica 
ceramic slip casting technique was used to form a 
ceramic shell of the models. The shell was then 
back-filled with a hydraulically setting magnesia 
ceramic for strength and support. Finally, the models 



Figure 2: MSL geometry. 


Parameter 

Full Scale 

Test Model 


1 wn ■ 

"" 1.56 

R„ (in.) 

38.79 

1.218 

1 R s (in.) 

3.89 

0.122 

R b (in.) 

13.95 

0.438 

L (in.) 

114.52 

3.596 

jj(deg) 

20 

20 

I (deg) 

70 

70 

C (deg) 

33.6 

33.6 


Table 1: Parameters for MSL geometry. 



Figure 3: MSL model and cavity orientation. 


were coated with a mixture of phosphors suspended in 
a silica-based colloidal binder. 

The proposed MSL entry vehicle is a 70-deg 
sphere-cone with a biconic afterbody. The geometry 
of the proposed MSL aeroshell is shown in Fig. 2 and 
its dimensions (full scale and test model) are shown in 
Table 1. The forebody will have six cavities of the 
same size and radial location spaced at 60-deg incre- 
ments, although the final size and locations have not 
been determined. 

The cast ceramic aeroheating models were 5-in. 
diameter, 0.0314-scale representations of a proposed 
13.29-ft diameter Mars Smart Lander aeroshell. Two 
model configurations were tested in this study. The 
first was a baseline configuration model without any 
cavities. The second configuration was a cavity para- 
metric model (Fig. 3). It had a total of six cavities 
located 60-degrees apart around the forebody. Three 
cavity diameters (w fs = 1.5-in, 2.2-in, and 3.0-in.) and 
two radial locations (r/R = 0.41 and 0.70) were 
included in this study. Each of the six cavities has 
been given a number and the radial locations and 
diameters of each are presented in Table 2. 


Cavity Number 

r/R 

w fs (in.) 

1 

0.70 

3.0 

2 

0.70 

2.2 

3 

0.70 

1.5 

4 

0.41 

3.0 

5 j 

0.41 

2.2 

6 

0.41 

1.5 


Table 2: Cavity numbering system for the MSL. 
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Test Matrix and Tunnel Conditions 


The data were collected on the baseline smooth 
model at angles-of-attack of 0-deg, 11-deg, 16-deg 
and 20-deg and freestream Reynolds numbers 
between 2.1xl0 6 and 5.8xl0 6 per foot. The cavity 
parametric models were tested at angles-of-attack of 
0-deg, 11 -deg, 16-deg and 20-deg at freestream Rey- 
nolds numbers between 2.4xl0 6 and 7.3xl0 6 per foot. 
The nominal test conditions are presented in Table 3. 

RESULTS AND DISCUSSION 

Flow Field Computations 



In order to refer to the cavities as model orienta- 
tion changes, two angles are defined. First, the model 
orientation angle, 0, is defined as 0-deg when cavity 1 
is on the vertical upright in an image, and increases as 
the model is rotated counterclockwise (Fig. 3; cavity 1 
is always on the +y' axis). The second angle, 6 , is the 
cavity orientation angle measured clockwise from the 
vertical upright in an image (6 = 0-deg on the vertical 
upright regardless of the value of 0). 

Small, circular marks were placed on the models 
to aid in data reduction and model orientation. These 
marks, referred to as fiducial marks, do not influence 
the flow over the model surface. The fiducial marks 
can be seen in run images as dark dots and should not 
be confused with cavities. 

Data Reduction 

One-dimensional, semi-infinite solid heat con- 
duction theory 11 was used to compute surface heating 
distributions from the global surface temperature data 
acquired through phosphor thermography. A constant 
heat-transfer coefficient is assumed in this theory, and 
empirical corrections 11 are made to account for 
changes in model substrate thermal properties with 
temperature. Phosphor images were acquired shortly 
after injection of the model to the tunnel centerline, 
which requires less than one second. 

Data cuts were extracted from the heat transfer 
images. Example data extraction locations are given 
in Fig. 4. Results are presented herein in terms of a 
non-dimensional heat transfer coefficient ratio, hlh FR , 
where h FR is the theoretical heating computed with 

the Fay-Riddell 13 method for a 1.2165-in. (3.09 cm) 
sphere, which is the radius of the spherical portion of 
the forebody of the test models, at 300 K. 

As detailed in Ref. 1 1 , the estimated uncertainty 
of the phosphor thermography technique is approxi- 
mately ±13% on the forebody. 


MSL flow field computations 4 were performed 
using the Langley Aerothermodynamic Upwind 

Relaxation Algorithm 5 ’ 6 (LAURA) at test facility con- 
ditions for laminar flow. A comparison between mea- 
sured experimental heating and computational heating 
predictions can be seen in Fig. 5. (A detailed discus- 
sion of boundary layer transition correlation and an 
aeroheating prediction analysis of the MSL for wind 
tunnel conditions is presented in the companion paper 
by Hollis 4 .) 

In order to determine MSL boundary layer prop- 
erties, an algorithm to was written to post-process the 
LAURA output files to extract boundary layer edge 


r 


1 


O 

J 

L. 


J 


a)ExperimentaI heating. 



b) Computational heating. 


0.0 0.1 0.2 0.3 




Figure 5: Comparison of experimental and com- 
putational results at Re^ = 3.0xl0 6 /ft, a = 16-deg. 
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Figure 6: Effect of angle-of-attack on boundary layer properties at Re „ = 3.0xl0 6 /ft. 


conditions. The boundary layer height (<5) was defined 
as the distance along the wall surface normal at which 
the stagnation enthalpy is equal to 99.5% of the 
freestream value. Flow properties were interpolated at 
discrete points along the wall surface normal and a 
numerical integration was performed to calculate the 
boundary layer momentum thickness (0). The stream- 
wise pressure gradient and Mach number were also 
calculated at the boundary layer edge. In the images 
presented in this section, all possible cavity locations 
are shown as white circles and streamlines are shown 
through all cavity locations. 

The effects of angle-of-attack are shown in Fig- 
ure 6 at Re oo = 3.0xl0 5 6 /ft. As angle-of-attack was 
increased, d and 0 on the windward surface of the 
forebody decreased, while on the leeward surface, 
they increased as expected. M e increased on the lee- 
ward surface of the forebody as angle-of-attack was 
increased and decreased on the windward surface, in 
general. The flow over the aeroshell forebody 
remained subsonic for the present range of 
angles-of-attack, and became supersonic only upon 
the flow expansion around the shoulder. The value of 
dp/ds increased along most of the symmetry plane 
with angle-of-attack (especially just downstream of 
the nose and near the stagnation point), but decreased 


just upstream of the nose and along rays approxi- 
mately 60-deg off of the windward centerline. 

The effects of freestream Reynolds number are 
shown in Figure 7 at a = 16-deg. As the freestream 
Reynolds number increased, 6 and 0 decreased over 
the entire forebody as expected. M e stayed approxi- 
mately constant with increasing freestream Reynolds 
number. The value of dp/ds decreased over the entire 
forebody except directly downstream of the nose as 
Re oo increased. 

Smooth Model Surface Heating Data 

The non-dimensional heating data for the base- 
line MSL model are shown along the model centerline 
in Figure 8 and the corresponding images are shown 
in Figure 9. The data are plotted versus the 
non-dimensional distance ratio r/R . The heating ratio, 
h/h FR , near the nose varied from approximately 0.68 
at a = 0-deg to approximately 0.72 at a = 20-deg. The 
heating ratio on the leeward surface of the forebody 
decreased to values ranging from 0.28 to 0.40 as 
angle-of-attack decreased from 20-deg to 0-deg and 
remained at these levels to the shoulder. Heating lev- 
els began to rise as the flow reached r/R of approxi- 
mately 0.9. The rapid expansion of the inviscid flow 
around the comer forms a large favorable pressure 
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gradient, which results in a reduction of the boundary 

layer thickness 14 . This increased heating near the 
shoulder increased with Reynolds number due to a 
further decrease in boundary layer height. 

The heating ratio on the conical portion of the 
windward surface decreased from the levels on the 
nose to approximately hfh FR = 0.50 for all non-zero 
angles-of-attack, then increased to between h/h FR = 
0.46 for a = 0-deg up to h/h FR = 0.75 for a = 20-deg 
near the shoulder due to expansion of the flow. 

Cavity Parametric Model Surface Heating Data 

The non-dimensional heating data extracted from 
the cavity parametric models were plotted beginning 
at the cavity centers and continuing along the center 
of the transitional/turbulent wedge which extended 
downstream from the cavity to the location of maxi- 
mum diameter (see Fig. 4; when a cavity did not pro- 
duce a wedge of increased heating, the data were 
taken along the cut location from a run that did trip 
the flow). The heating data are plotted verses the 
non-dimensional distance s/R , which is the distance 
from the center of the cavity over the surface of the 
model to the outer diameter of the model non-dimen- 
sionalized by the model radius. Included with every 
cavity plot is a comparison with data extracted along 
the same line from the baseline smooth model at Re & 

= 3.0xl0 6 /ft. 

The presence of a cavity may result in high, 
localized heating at the downstream edge of the cavity 
due to flow separation and reattachment within the 
cavity, and may be accompanied by viscous effects 
downstream of its location (i.e. boundary layer transi- 
tion). A cavity’s ability to influence these parameters 
will be referred to as its effectiveness. 

The effect of cavity diameter at three Reynolds 
numbers is shown in Fig. 10. As expected, the effec- 
tiveness of the cavities is shown to increase as Re & 
increases. This is because the height of the boundary 
layer decreased as Re x increased (Fig. 7). Also as 
expected, the larger the cavity size, the greater the 
effect on the state of the boundary layer. The effec- 
tiveness of the cavities ranged from essentially zero at 
the lower values of Re & to an increase in downstream 
heating levels of approximately 325% for Re ^ = 

7.3xl0 6 /ft (see Ref. 4 for comparison to computed 
turbulent values). 

As angle-of-attack was increased (Fig. 11), the 
effectiveness of a cavity was strongly influenced by its 
location on the forebody. Cavities at both radial loca- 
tions for 6 =1 20-deg and 240-deg became more effec- 
tive with increasing angle-of-attack; and also at r/R = 


= 0-deg 


Approximate data 
extraction line 



0 = 240-deg 


= 60-deg 


= 1 20-deg 


0= 180-deg 



0.40 0.60 

s/R 

(C) R9 x = 7.3x1(fm 


Figure 10: Effect of cavity diameter at a = 16-deg 
and 0 = 240-deg for Re„ = 3.4xl0 6 , 5.1xl0 6 , and 
73xl0 6 /ft. 
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low (compared to the second pattern) and the down- 
stream pressure gradients were strongly negative. The 
second trend (Fig. 12.f) was an immediate, relatively 
small increase in heating downstream of the cavity 
followed by a gradual increase in heating levels. Cavi- 
ties exhibiting this behavior were located at both 
radial locations on the leeward centerline and at r/R = 
0.70 at 6= 60-deg and 300-deg at all angles-of- attack. 
The r/R = 0.41, 6 = 60-deg and 300-deg locations 
were included in this pattern at a = 11 -deg. These 
locations correspond to regions where 6 and 0 were 
relatively large (compared to the first pattern) and the 
downstream pressure gradients were greater than 
those of the first pattern. 

The present measurements in the 20-Inch Mach 6 
Air Tunnel for the proposed MSL aeroshell demon- 
strated that the effectiveness of a given cavity depends 
on the local values of d and 0. As these parameters 
decreased, the effectiveness of the cavity was shown 
to increase. One method of predicting boundary layer 
transition due to the presence of a cavity is to compare 
the Reynolds number based on momentum thickness 
with the ratio of cavity diameter to boundary layer 
height (w IS). One can expect the effectiveness of a 
constant cavity size to increase with w/<5, as shown. 
This correlation is discussed in Ref. 4 for the pro- 
posed MSL aeroshell. 

CONCLUSIONS 

An experimental aeroheating study was con- 
ducted on a proposed Mars Smart Lander aeroshell in 
the NASA Langley 20-Inch Mach 6 Air Tunnel using 
the technique of phosphor thermography. This test 
was conducted in order to determine if heat shield 
cavities of significant size have an impact on forebody 
heating. The study concluded that forebody cavities 
did indeed have a profound impact on heating levels 
and should be considered in the design process of the 
Mars Smart Lander. 

This study was complemented by an analysis of 
boundary layer quantities obtained from thin-layer, 
laminar, Navier-Stokes flowfield solutions. The com- 
putational results quantified the boundary layer height 
and momentum thickness as functions of freestream 
Reynolds number, angle-of-attack, and location on the 
aeroshell forebody. Tests were conducted on baseline 
smooth and cavity parametric Mars Smart Lander 
aeroshell models at freestream Reynolds numbers 
from 2.6x1 0 6 /ft to 7.3xl0 6 /ft. 

Smooth model heating data showed that heating 
levels near the shoulder increased due to the thinning 
of the boundary layer. This increase in heating near 


the shoulder increased with Reynolds number. At an 
angle-of-attack of 20-deg at a Reynolds number of 

5.8xl0 6 /ft, natural transition raised heating levels on 
the windward surface downstream of the stagnation 
point by approximately 14% on the windward surface. 

Cavity parametric model heating data showed 
that the effectiveness of a cavity depended on the cav- 
ity diameter and boundary layer height. As the cavity 
diameter increased, the effectiveness of the cavity 
increased. As boundary layer height decreased, the 
effectiveness of a cavity increased. The presence of a 
cavity was shown to cause the increase in downstream 
heating levels by as much as 325% due to boundary 
layer transition. 

The present study has provided a comprehensive 
data base (approximately 800 data points) on the state 
of the boundary layer of a proposed MSL aeroshell 
with cavities located on the forebody heat shield. This 
database was used to formulate a transition correla- 
tion to predict the onset of boundary layer transition 
in a companion paper by Hollis 4 . 
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Re*, 

(1/ft) 

M x 

T x 

(K) 

Poo 

(kg/m 3 ) 

Voo 

(m/s) 

Rl’oc 

(1/m) 

h F R 

(kg/m 2 -s) 

Ifr 

(W/cm 2 ) 

2.1xl0 6 

5.95 

62.0 

3. 35x1 0' 2 

938x6 

7.07x1 0 6 

0.283 

5.70 

2.6xl0 6 

5.97 

62.2 

4.05xl0' 2 

943.0 

8.55xl0 6 

0.313 

6.45 

3.0xl0 6 

5.98 

62.2 

4.62x1 0" 2 

944.4 

9.79xl0 6 

0.335 

6.95 

3.4xl0 6 

5.99 

61.6 

5. 29x1 O' 2 

940.1 

1.13xl0 7 

0.356 

7.25 

42xl0 6 

6.00 

61.3 

6.41xl0' 2 

940.4 

1.37xl0 7 

0.392 

7.95 

4.6xl0 6 

6.01 

63.5 

7. 25x1 O' 2 

958.5 

1.52xl0 7 

0.427 

9.48 

5.1xl0 6 

6.02 

63.4 

7.92x1 O' 2 

958.7 

1.67xl0 7 

0.446 

9.92 

5.8xl0 6 

6.03 

62.8 

8. 99x1 O' 2 

955.6 

1.91xl0 7 

0.474 

10.4 

73xl0 6 

6.06 

62.3 

1. 13x10'* 

954.6 

2.41xl0 7 

0.529 

11.5 


Table 3: Nominal flow conditions of the Langley 20-Inch Mach 6 Air Tunnel. 
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